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Abstract 

Given two polynomials, we find a convergence property of the CCD of the rising factorial 
and the falling factorial. Based on this property, we present a unified approach to computing 
the universal denominators as given by Gosper's algorithm and Abramov's algorithm for 
finding rational solutions to linear difference equations with polynomial coefficients. 
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1 Introduction 

Let N be the set of nonnegative integers, K be a field of characteristic zero, ]K(n) be the 
field of rational functions over IK, and ]K[n] be the ring of polynomials over K. We assume 
that subject to normalization the gcd (greatest common divisor) of two polynomials always 
takes a value as a monic polynomial, namely, polynomials with the leading coefficient being 
1. Recall that a nonzero term t„ is called a hypergeometric term over K if there exists a 
rational function r € K(n) such that 



t 



n+l 



r{n). 



If r(n) = a{n)/b{n), where a{n), h{n) G K[n], then the function a{n)/b{n) is called a rational 
representation of the rational function r(n). If gcd(a(n), = 1 holds, then a{n)/b{n) is 
called a reduced rational representation of r(n). 

Gosper's algorithm [6] (also see [TJ [HI [HI [T71 [18l [191 [16] ) has been extensively studied and 
widely used to prove hypergeometric identities. Given a hypergeometric term Gosper's 
algorithm is a procedure to find a hypergeometric term satisfying 

Zn+l — Zn = tm (1-1) 



if it exists, or confirm the nonexistence of any solution of (II. ID . The key idea of Gosper's 
algorithm lies in a representation of rational functions called Gosper representation; i.e.. 



writing the rational function r(n) in the fohowing form: 

, , a(n) c(n + 1) 
b[n) c[n) 

where a, b and c are polynomials over IK and 

gcd(a(n), 6(n + h)) = 1 for all /i G N. 

Petkovsek [12] has realized that a Gosper representation becomes unique, which is called the 
Gosper-Petkovsek representation, or GP representation, for short, if we further require that 
b, c are monic polynomials such that 

gcd(a(n), c(n)) = gcd(6(n), c(n + 1)) = 1. 



In the same paper, Petkovsek also gave an algorithm to compute GP representations; 
subsequently we will call it the "GP algorithm" . In [Tl] , Paule and Strehl gave a derivation 
of Gosper's algorithm by using the GP representation. In [10], equipped with the Greatest 
Factorial Factorization (GFF), Paule presented a new approach to indefinite hypergeometric 
summation which leads to the same algorithm as Gosper's, but in a new setting. In [9], 
Lisonek and et al, gave a detailed study of the degree setting for Gosper's algorithm. 

Finding rational solutions is important in computer algebra because many problems can 
be reduced to rational solutions. For example, we may consider the generalization of Gosper's 
algorithm. Given a linear difference equation 

d 

Pm{n)y{n + m) = p(n), (1.2) 

m=0 

where po(^)) Pi • • • TPd{n),p{n) € K[n] are given polynomials such that po{n) ^ 0, Pd{n) ^ 
0, a polynomial g{n) G ]K[n] is called a universal denominator for (jl.2p if and only if for 
every solution y{n) € K.(n) to ()1.2p there exists a f{n) G ]K[n] such that y{n) = f{n)/g{n). 
Once a universal denominator is found, then it is easy to find the rational solutions of the 
linear difference equation (|1.2p by finding the polynomial solutions using the techniques in 
[21 Ul |T2] . Abramov [2] developed an algorithm to find a universal denominator of (jl.2p which 
relies on all the coefficients Po{n),pi{n), . . . ,pd{n),p{n). In [3], an improved version is given 
which requires only two coefficients po{n) and Pd{n). Compared with the simplicity of the 
output of Abramov's algorithm, the justification is quite involved. Recall that the dispersion 
dis(a(n), 6(n)) of the polynomials a(n),6(n) G K[n] is the greatest nonnegative integer k (if 
it exists) such that a(n) and b{n + k) have a nontrivial common divisor, i.e., 

dis(a, b) = max{/c G N | deg gcd(a(n), b{n + A;)) > 1}. 

If k does not exist then we set dis(a,6) = —1. Observe that dis(a(n), 6(n)) can be com- 
puted as the largest nonnegative integer root of the polynomial R{h) G M.[h] where R{h) = 
Resn{a{n), b{n + h)). 

The main result of this paper is the discovery of a convergence property of the GCD of 
rising factorial of a polynomial b{n) and the falling factorial of another polynomial a(n). By 
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using the limit of the GCD sequence, we may transform a rational difference equation into 
a polynomial difference equation. The convergence argument yields a new and streamlined 
approach to the explicit formula for Abramov's universal denominator. Note that this explicit 
formula can be used to compute rational solutions of a linear difference equation ()1.2p . In 
addition, we derive Abramov's universal denominator from Barkatou's explicit formula. The 
relation between Barkatou's approach and Abramov's algorithm has been discussed in detail 
by Weixlbaumer [15j . 

2 The Convergence Property 

The main idea of this paper is the following convergence property of a GCD sequence. It 
turns out that this simple observation plays a fundamental role in finding rational solutions of 
linear difference equations, and it can be viewed as a unified approach to several well-known 
algorithms. 

Theorem 2.1. Let a{n) and h{n) he two nonzero polynomials in n and let 

ko = dis(a(n - 1), 6(n)) = max{A; G N | deg gcd(a(n - 1), b{n + k)) > 1}. (2.1) 

Define 

Gkin) = gcd(6(n)6(n + I) . . .b{n + k - l),a{n - l)a{n - 2) . . . a{n - k)). (2.2) 
Then the sequence Gi(n),G2(n), . . . converges to Gko+iin). 

Proof. For all k > ko we have 

gcd(a(n -l),b{n + k)) = ■■■ = gcd(a(n -k-l),b{n + k)) = 1. (2.3) 

Note that 

Gk+i{n) = gcd(6(n)6(n + 1) . . . 6(n + k),a{n - l)a(n -2)...a{n-k- 1)). 
This implies that 

Gfc(n) = Gfc+i(n), 

for all k > kf). Moreover, one sees that once (j2.3p is satisfied for > /cq, it is also satisfied for 
k + 1. It follows that 

Gfc(n) = Gfc+i(n) = Gk+2{n) = ■■■ , 
for all k > kf), and this completes the proof. □ 

Using the above convergence property of the sequence Gkin), we are led to a simple 
approach to Gosper's algorithm without resorting to the Gosper representation or GP rep- 
resentation of rational functions. Given a hypergeometric term t„ and suppose that there 
exists a hypergeometric term Zn satisfying equation (II. ip . then by using (II. ip we find 

r(n)2/(n + 1) - y(n) = 1, (2.4) 

where r(n) = tn+i/tn and y{n) = Zn/tn are rational functions of n, see [H]. 
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Theorem 2.2. Let r{n) and y{n) in equation l\2.4\ l be in terms of their reduced rational 
representations: 

h[n) g(n) 

Then 

g{n) I Gfeo+i(n), 
where and Gk{n) are defined in Theorem \2.1[ 



Proof. Using (12. Sp in (12. 4p gives 

a{n)g{n)f{n + 1) - b{n)g{n + l)/(n) = b{n)g{n)g{n + 1). (2.6) 
From the above relation, we immediately get that 

g{n) I b{n)g{n + 1) and that g{n + 1) | a{n)g{n). 
Using these two relations repeatedly we obtain 

g{n) \ b[n)b[n + 1) . . . b{n + k — l)g{n + k), 
g{n) j a(n — l)a(n — 2) . . . a{n — k)g{n — k), 
for all k ^'H. Since IK has characteristic zero, 

gc<l{g{n),g{n + k)) = gcd{g{n),g{n - k)) = 1, 
for all large enough k. It follows that 

g{n) I b{n)b{n + 1) . . . 6(n + /c - 1), (2.7) 
g{n) I a{n - l)a(n - 2) . . . a(n - k), (2.8) 
for all large enough k. Therefore 

g{n) I Gfc(n), 

for all large enough k. The rest of the proof follows when k goes to infinity in this equation 
and by Theorem 12.11 □ 

The next step is simply to set 

g{n) = Gk^+i{n) (2.9) 

in equation (j2.6p as in the GFF algorithm of Paule. If equation (j2.6p can be solved for 
f eK[n], then 

is a hypergeometric solution of (jl.ip : Otherwise no hypergeometric solution of (jl.ip exists. 
Note that the solution /(n) may not be coprime to g{n). However, it is clear that this does 
not affect the solution of y(n). Indeed, the polynomials /(n) and g{n) can be recovered from 
the solution of y{n) after dividing the greatest common factors. 
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Algorithm 2.3. 

INPUT: r{n) G ]K(n) such that tn+i/tn = r{n) for large enough n in N. 

OUTPUT: a hypergeometric solution Zn of jp if it exists, otherwise "no hypergeometric 
solution of (jl.ip exists". 

(1) Decompose r{n) into a/b where a, b are two relatively prime polynomials. 

(2) Compute as in <\2.1\ . 

(3) If > then compute g{n) = Gko+iin), where Gk{n) is defined as in l\2.^ . otherwise 
gin) = 1. 

(4) // equation can he solved for f G ]fC[n] then return Zn = i™; Otherwise return 
"no hypergeometric solution of (|i.ip exists". 

Let us take an example from |14j : 
Example 2.4. Let t^ = (4n + 1) • (2n+i)\ ' 

4n + 5 

^ ' tn 2(4n + l)(2n + 3)' 

Hence a{n) = 4n + 5, 6(n) = 2(4n + l)(2n + 3) and then k^ = 0. Note that for all k > fco, 
equation (l^.gp zs satisfied. From (l^.gp . = n + |. -By (j^.gp . /(n) is a polynomial which 
satisfies 

2/(n + 1) - 4(2n + 3)/(n) = (2n + 3)(4n + 1). 
The polynomial f{n) = —^[2n + 1) is a solution of this equation. Therefore, 

f{n) n! 

We remark that the argument for the relations (|2.7p and (|2.8p is used by Petkovsek |12j . 
Moreover, the products on the right hand sides of ()2.7p and (j2.8p can be written in the 
notation of rising or falling factorials as introduced by Paule [TO] , 



3 Connections to Gosper's and Abramov's Algorithms 

We will show how Theorem 12.21 is related to Gosper's algorithm and Abramov's algorithm 
for finding rational solutions of linear difference equations with polynomial coefficients [3] . 

Abramov's Algorithm (general order d): Consider the difference equation 

Pd{n)y{n + d) + . . . + po{n)y{n) = p{n) (3.1) 

with given po{n),pi{n), . . . ,P(i{n),p{n) S K[n] such that po and pa are nonzero. Abramov 
gave the following algorithm to compute a universal denominator G{n) for (j3.ip : Define 

= dis{pd{n — d),po{n)) = max{A; € N | deg gcd{pd{n — d),pQ{n + k)) > 1}. (3.2) 
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If = 
define 



1 set G(n) = 1, i.e., in this case all rational solutions are polynomials. If > 0, 

(3.3) 



AN+i{n) = pd{n - d), BN+i{n) = po{n), 
and for i = N down to i = do: 



di{n) 
Ai{n) 



gcd(^i+i(n),5i+i(n + i)), 
Ai+i{n) , „ , . Bi+i{n) 



and Bj (n) 



di{n) di{n-i)' 
If we use the notation of the falling factorial of a polynomial introduced in [10] by 

[f[x)f=f{x)f{x-l)---f{x-k + l), 
then Abramov's universal denominator of (j3.ip can be written as 

G{n) = [do(n)]i[di(n)]^. • • [d^in)]^. 
There is also an explicit formula for Abramov's universal denominator ()3.6p . namely, 

G{n) = gcd([po(n + N)]^, [p,(n - d)]^). 



(3.4) 
(3.5) 



(3.6) 
(3.7) 



It can be seen that (13. 6p is equivalent to Theorem 3 in Abramov-Petkovsek-Ryabenko [ij. A 
generalized form of (13. 7p can be found in Barkatou [5]. For completeness, we give a proof of 
the fact that the presentations ()3.6p and (j3.7p indeed coincide with each other. To this end 
we will follow the survey of Weixlbaumer [15] . 



First of all, based on the definition of A^ and the fact that Ai \ ^j+i and Bi \ Si+i, it 
follows that for 0<i<A^ + l, we have 



gcd{Ai{n),Bi{n + k)) = I for all k> N. 
Moreover, from ()3.4p and (j3.5p . for < i < A^ it follows that 

gcd{Aj{n),Bj{n + i)) = I for < j < i. 



(3.8) 
(3.9) 



Therefore, 



G{n) 



gcd(bo(n + A^)]^, [Pd{n - d)]^) 
god{[BM+i{n + AT)]^, [AN+i{n)]^) 



gcd 



[dN{n)] 



BN+i{n + N)- 


N+l 


' AN^i{n)' 


N+l\ 


dN{n) 




dN{n) 





[dN{n)] 



N+l 



N+l 



.cd{[BM{n + N)]^, [AN{n)] 



Observe that 



ycd{BN{n + N),[AN{n)]^) 

= gcd{BN{n + N),AN{n)) 



[by ^ 

[by m 



and similarly, 

gcd([S^(n + N- 1)]^, AN{n - N)) = 1. 

Consequently, 

gcd([B^(n + N)]^, [AN{n)]^) = gcd{[Br,{n + N- 1)]^, [^7v(n)]^). 
In the same manner one can successively split off the factors 

[dw_i(n)F,...,[do(n)]i 
until one arrives at (j3.6p , which completes the proof of the equality of (j3.6p and (j3.7p . 



Next we remark that Theorem 12.21 follows from Abramov's algorithm, strictly speaking, 
the universal denominator given by Abbramov's algorithm. Equation (|2.4p is equivalent to 



a{n)y{n + 1) — b{n)y{n) = b{n), (3.10) 

which is (|3.ip with d = 1, pi{n) = a(n), po{n) = —b{n), and p{n) = h{n). From ()3.7p . 
Abramov's algorithm gives the following universal denominator of (j3.10p 

G{n) = gcd([a(n - 1)]^, [b{n + iV)]^). (3.11) 



where N = hy (j2.ip . Using (j2.2p we have that G{n) = GfcQ+i(n), hence Theorem [27 
determines the same universal denominator as Abramov's algorithm. 

Next we show that Abramov's algorithm delivers a Gosper representation for r(n) = 
tn+i/tn if r{n) = a[n)/h{n) is the reduced rational representation of r(n). Prom ()3.5p we 
obtain that a{n — 1) = ^7v+i("-) = dN{n)Ai~j{n), and by iteration, 

a{n) = do{n + 1) di(n + 1) . . . dN{n + 1) Ao{n + 1). (3.12) 

Analogously, ()3.5p implies that 

6(n) = -do{n) di(n - 1) . . . (iiv(n - N) Bo{n). (3.13) 

Consequently, in view of representation (j3.6p for G{n), one obtains that 

a(n) ^ G(n + 1) Ao(n + l) 

6(n) G(n) Bo{n) ' ^ ^ 

Note that in (j3.10p w.l.o.g. we can assume that gcd(a(n), 6(n)) = 1, which then implies 

gcd{Aoin + l),Boin)) = l. (3.15) 
But more is true. Namely, by (j3.9p 

gcd(Ao(n + 1), Bo{n + i + 1)) = 1 for < i < iV, 

i.e., 

gcd{Ao{n + l),Bo{n + i)) = 1 forl<i<A^ + l, (3.16) 



and by (jHTSD . 

gcd(Ao(n + l),Bo{n + A; + 1)) = 1 for ah k> N, 



I.e., 

gcd{Ao{n + l),Bo{n + k)) = 1 for all A: > iV + 1. (3.17) 

Finally, combining (j3.15p . (j3.16p and (|3.17p into one condition results in 

gcd(Ao(n + l),So(n + /i)) = 1 for all /i > 0. (3.18) 
Hence the right hand side of (j3.14p is a Gosper representation for a(n)/6(n). 
Example 3.1. Let = ^"^^^ , then 

tn+i a{n) 



r{n) 



6(n)' 



where a{n) = n + 2>, h{n) = [n + l)(n + 2). From Ahramov's algorithm, we have N = 1. 
By using Abramov's algorithm with AN+i{n) = A2{n) = a{n — 1) = re + 2 and BN^i{n) = 
B2{n) = h{n) = (n + l)(re + 2), we obtain that AqIji) = 1, Bo(n) = n + 2, and the universal 
denominator G{n) = (re + l)(n + 2). Note that 

gcd(G(re + l),5o(re)) =re + 2. 

This means that the Gosper representation ^3A^ in general is not the GP representation for 
a{n)/b{n). 



As explained in Paule [10], also the GP algorithm [12] for finding a GP representation of 
r(n) = a{n)/b{n) computes a universal denominator for ()3.10p . namely as follows. 

Petkovsek's GP Algorithm: Compute N as in Abramov's algorithm; i.e., N = Uq as, in 
(I3TT]) . liN = -l set n(re) = 1. If iV > define 

ao(re) = a(re), 6o("-) = ^in), (3.19) 

and for i = luptoi = A^+l do: 

5i{n) = gcd(ai_i(n),6i_i(re + i)), (3.20) 

•^A^j = . / X and bi{n) = — -. (3.21) 

di[n) di[n-z) 

This determines a universal denominator by setting 

u{n) = Mn - l)]HHn - 1)]^- • • [SN+i{n - 1)]^ . (3.22) 

Note that this algorithm essentially consists in running the loop in Abramov's algorithm in 
the REVERSE direction. (Note that also its initialization is slightly different, namely starting 
with the pair (a(re),6(re)) instead of (a(n — l),6(n))). 



Finally, as above, from ()3.2ip and by using ()3.22p one obtains that 

a(re) u(n + 1) aAr+i(re) 
b{n) u{n) bN+i{n)' 



(3.23) 
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Petkovsek's algorithm is designed in such a way that (j3.23p is not only a Gosper represen- 
tation, but also a GP representation for a{n) /h{n); see |12j . |14j . and also [10]. This means 
that besides 

gcd(aAr+i(n), 6iv+i(n + h)) = 1, for all /i > 0. (3.24) 

we also have that 

gcd(u(n + l),6^+i(n)) = 1, (3.25) 

and 

gcd(ti(n),aAr+i(n)) = 1. (3.26) 

From ([XTill and (13:23]) we get 



(3.27) 



u{n + l) aN+i{n) _ G{n + 1) Ao(n + 1) 
u{n) bN+i{n) G{n) Bo{n) 

By Lemma 5.3.1 (see |14J, P-82), we obtain that 

u{n) I G{n), 

which implies that the universal denominator given by GP algorithm is a factor of the uni- 
versal denominator given by Abramov's algorithm. 



4 Rational Solutions of Linear Difference Equations 

In this section we show how to deduce the explicit formula (|3.7p for Abramov's universal 
denominator by using the convergence argument. 

Theorem 4.1. Let poin) and Pd{n) he two nonzero polynomials in n, and let N he defined 
as in ^SIM- Put 

(k-l k-l \ 

llpo{n+j),Ylpd{n-d-j)\. (4.1) 
i=o j=o J 

Then the sequence Gi{n),G2{n), . . . converges to GN+i{n). 

Proof. For all A; > we have 

gcd{po{n),pd{n -d-k)) = ■■■ = gcd(po(n + k),pd{n - d - k)) = I. (4.2) 

Note that 

(k k \ 

J|po(ra + j), JJpd(n -d- j) \ . 
3=0 j=o J 

This implies that 

Gfc(n) = Gfc+i(n), 

for all k > N. Moreover, one sees that once (|4.2p is satisfied for k > N, it is also satisfied for 
k + 1. It follows that 

Gfc(n) = Gk+i{n) = Gk+2{n) = ■■■ , 
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for all A; > A^, and this completes the proof. □ 

By using the convergence property of the sequence Gfc(n), we can obtain the explicit 
formula (|3.7p for Abramov's universal denominator. Our proof requires only 1cm and gcd 
computations. 

Theorem 4.2. Given a linear difference equation 

d 

'^Pm{n)y{n + m) = p{n), (4.3) 

m=0 

where po{n),pi{n), . . . ,pd{n),p{n) G K[n] are given polynomials such that po^n) ^ 0, pd{n) ^ 
0. Let f{n)/g{n) he the reduced rational representation of y{n), and N he the dispersion of 
Pd{n — d) andpo{n). Then 

G{n) = gcd([po(n + iV)]^, [p,(n - d)]^), 

as given hy 7[ ), is a universal denominator of rational solutions of 

Proof. From (14. Sp it follows that 



E j[n-\-m 
Vm{n) ; ( =p{n). 4.4 
qin + m) 
m=o ' 



f{n + m 

PmKn; 

m=0 

Letting 

l{n) = lcm{g{n + l),g{n + 2),...,g{n + d)), (4.5) 
and multiplying equation ()4.4p by l{n)g{n), we obtain 

d ^\ 

Pm{n)f{n + m)— — ■ g{n) = p{n)l{n)g{n). (4.6) 

^-^ qin + m) 

m=0 UK ' I 

From (j4.5p . we have the following divisibility conditions: 

g{n + m) | /(n) for m = 1, 2, . . . , d. 

Thus l{n) / g{n + m) are polynomials for m = 1, 2, . . . , d. From (|4.6p we obtain 

5(n) |po(f^) •lcm(5r(n + l),5r(n + 2),...,5r(n + d)). (4.7) 

Similarly, multiplying equation (j4.4p by /(n — l)5'(ra + d) and then substituting n — d for n, 
we obtain that 

g{n) I pd(n - d) • lcm.{g{n - l),g{n -2),...,g{n- d)). (4.8) 
Shifting n by 1 in ()4.7p yields 

5(n + l) |po(f^+l) •lcm(5(n + 2),5(n + 3),...,(7(n + d+l)). (4.9) 
Substituting (j4.9p into (j4.7p we see that g{n) divides 

Po(n) ■ lcm{po{n + 1) • lcm(5f(n + 2), ^(n + 3), . . . , 5r(n + d + l)),g{n + 2),...,g{n + d)). 
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So we can write 

g{n) j po{n)po{n + 1) • lcm{g{n + 2), g{n + 3), . . . , g{n + d + 1)). 
By induction we may derive for k > 1, 

5((n)| JJpo(f^ + j) • lcin{g{n + k),g{n + k + 1), . . . , g{n + k + d - 1)). 

j=0 

It foUows that 

g{n)^ n ?'o("' + •?') ■ + ^)5'("' + + 1) . . . + k + d - I). 

j=0 

Since IK has characteristic zero, there is a large enough k such that for any j > k 

gcd{g{n),g{n + j)) = 1. 

It fohows that 

k-l 

5-(n)| J|po(f^ + j), 

j=0 

for ah large enough k. Analogously, from ()4.8p we get 



k-l 

c/(n)| Y[pd{n-d-j), 

j=0 

for all large enough k. Therefore 

g{n) I Gfe(n), 

for all large enough k, where Gkin) is defined as in (j4.ip . Setting A; to infinity in this equation, 
by Theorem 14. II we get 

(N N \ 

'[{po{n + j),Y{Pdin-d-j)\ =G{n), 
j=o j=o J 

as desired. □ 



From equation (j4.4p we get 

d d d 

^Vm{n) ■ f{n + m) ■ g{n + j) = p{n) ■ Y\_9{n + j)- (4.10) 

m=0 

The next step is simply to set 



m=0 j=o j=o 



g{n) = G{n), 

in equation (I4.10|) . If equation ()4.10p can be solved for / G K[n] then y{n) = f{n)/g{n) is a 
solution of (14.31): otherwise (14.31) has no rational solutions. 
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Algorithm 4.3. 

INPUT: nonzero polynomials Po{n),pd{n). 
OUTPUT: a universal denominator g{n) of 



(1) Compute N = dis(pd(n — d),pQ{n)) = max{k € N | deg gcd(pd(n — d),pQ{n + k)) > 1}. 

(2) If N > then compute g{n) = G{n), where G{n) is defined as in l \3.7\) , otherwise 
g{n) = 1. 

Example 4.4. Find a rational solution of the equation 

(n + 4)(2n + l)(n + 2)y{n + 3) - (2n + 3)(n + 3)(n + l)y{n + 2) 

+n(n + 2)(2n - 3)y{n + 1) - (n - l)(2n - l)(n + l)y(n) = 0. (4.11) 

M^e have po{n) = -(n - l)(2n - l)(n + 1), prf(n) = (n + 4)(2n + l)(n + 2), i/ien iV = 2 and 
i/ien g{n) = (n — l)(n + l)n. By l\4.10\) , fin) is a polynomial which satisfies 

n{2n + l)(n + 2)(n + l)/(n + 3) - n(2n + 3)(n + 3)(n + l)/(n + 2) + n{2n - 3) 
•(n + 3)(n + 2)/(n + 1) - (2n - l)(n + 3)(n + l)(n + 2)/(n) = 0. 

T/ie polynomial f{n) = Cn{2n — 3) is a solution of this equation. Thus 

f{n) (2n + l) 
y{n) = — — = C- 



c/(n) (n2 - 1) 
is a rational solution of 
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